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Periodic boundary conditions have not a unique implementation in magnetic systems where all 
spins interact with each other through a power law decaying interaction of the form 1/r", r being 
the distance between spins. In this work we present a comparative study of the finite size effects 
oberved in numerical simulations by using first image convention and full infinite of periodic bound- 
ary conditions in one and two-dimensional spin systems with those type of interactions, including 
the ferromagnetic, antiferromagnetic and competitive interactions cases. Our results show no sig- 
nificative differences between the finite size effects produced by both types of boundary conditions 
when the low temperature phase has zero global magnetization, while it depends on the ratio ct/d 
for systems with a low temperature ferromagnetic phase. In the last case the first image convention 
gives much more stronger finite size effects than the other when the system enters into the classical 
regime a/d < 3/2. 
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I. INTRODUCTION 

Boundary conditions are a central issue in almost every statistical mechanics calculation. In particular, periodic 
boundary conditions are the usual way of diminishing finite size border effects in translationally invariant systems. 
This point is of special importance in numerical simulations, where an appropriated handhng of finite size effects 
may determine completely the quality of the results. While both the interpretation and implementation of periodic 
boundary conditions are straigthforward in systems with short range microscopic interactions, their usage in systems 
with long-range interactions is not so easy, with several subtleties to take into account. In this work we will consider 
magnetic systems where the iteractions between pairs of spins decay as l/r", r being the distance between spins and 
a > (this includes, among the most conspicuous examples. Coulomb and dipolar interactions). Then, every spin 
interacts with any other spin in the system, with a slow decaying distance-dependent intensity. Finite size effects 
are much more stronger in these kind of systems than in systems with short range interactions, because every spin 
feels directly the infiuence of the border. Hence, boundary conditions are a topic of central importance in numerical 
simulations of these kind of systems. 

Spin systems with the above type of interactions may present a very rich phenomenology, even in one dimension, 
and have attracted a lot of attention in the last years, both because its applications to ultrathin magnetic films 
systems 01, glassy systems 0| and to the study of statistical mechanics fundamentals among many others. 

There are two alternative representations of a finite system with periodic boundary conditions. Let us think on 
a spin system defined on a hypercubic d-dimensional lattice with N = L'^ sites. On one hand, we can visuaHze the 
system as a d-torus (closed ring in d = 1; torus in d = 2, etc.), where a given spin interacts with its closer neighbors 
defined by the topology of the d-torus, up to a certain range of interaction. On the other hand, we can think that 
we have an infinite system, where the original finite system has been replicated infinite times in all the coordinate 
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directions. In other words, we can think that the infinite lattice has been partitioned into cells of size N = L"^ and 
that we have chosen only the periodic solutions of the problem, with periodicity L in all the coordinate directions. 

While the difference between both views is just a matter of interpretation for systems with short-range interactions, 
the situation changes for long-range interacting systems. According to the second scheme, a given spin will interact 
with infinite replicas of everyone of the rest of the spins and we can express the effective interaction between two 
spins inside the system as an infinite sum over replicas (this includes a "self-interaction" term, that accounts for the 
interaction of every spin with its own replicas). We will call these infinite periodic boundary conditions (IPBC). On 
the other hand, in the closed topology of the first scheme it is consistent to consider that every pair of spins inside the 
system interacts only through its minimal distance over the hypertorus. This corresponds to the minimal truncation 
of the infinite series of the previous scheme and it is sometimes called the first image convention; we will refer to it 
as first image periodic boundary conditions (FIPBC). While both type of boundary conditions give the same result in 
the thermodynamic limit L ^ oo, they are clearly different for finite systems. 

Now, the question is how good are both type of boundary conditions to wipe out finite size effects in systems of 
moderate sizes, such as those available in practical numerical calculations? This is of particular interest, considering 
that the numerical simulations are much more costly than in the case of short range interactions (the computational 
complexity of the algorithms is typically 0(N'^)). Although both IPBC and FIPBC have been used in several works, 
including both numerical 0-001 and analytical^3^3| calculations, to the best of our knowledge the only comparison 
of its relative efficiency in magnetic systems was reported in an old work by Kretschmer and Binder [iJl, in the case 
of dipolar interactions in a three-dimensional system. In that work the authors observed a very poor performance 
of the FIPBC compared to the IPBC when the system exhibits a ferromagnetic phase. Since then, several authors 
assumed this result as a general rule, independent of the dimension d, the range of the interactions (i.e., of the 
exponent a) and the symmetry of the low temperature phase. In this work we compare the relative efficiency of both 
type of boundary conditions by considering the equilibrium critical behaviour of ferromagnetic, antiferromagnetic and 
competing interactions in d = 1 and d = 2 and for different values of the exponent a, using a single spin-fiip Monte 
Carlo algorithm (Metropolis). 

While the implementation of the FIPBC is straightforward, the implementation of the IPBC is much more cum- 
bersome and needs care, because the infinite series for the effective interactions are usually slowly convergent. This 
is another reason for studying its relative efficiency. The usual way of handling numerically the IPBC in d = 2 and 
d = 3 is to adapt the Ewald sums technique, originally derived for systems of interacting charged na,rticles[Tfll \2(\. 
to magnetic models[l, 18). In section^ we compare the critical temperature of one dimensional Ising models with 
ferromagnetic, antiferromagnetic and competing interactions, for different system sizes and 1 < a < 2 using FIPBC 
and IPBC. In this case we introduce a numerical technique for implementing the IPBC that is easier to implement 
in d = 1 for arbitrary values of a than the Ewald sums technique. We also analized the critical behaviour of the 
antiferromagnetic model for < a < 1. In section UTTI we compare the specific heat numerical results obtained around 
the critical region for ferromagnetic long range interactions with a = 3, and competing ferromagnetic short range 
interactions with antiferromagnetic (dipolar) long range interactions, using the Ewald sums for the IPBC. A general 
discussion of the results is presented in section Hvl 

II. ONE-DIMENSIONAL ISING MODELS 

We considered first an Ising model of the type 

H = -JoY, finj) S,S, = ±1 (1) 

where the sum j) runs over all distinct pairs of sites (i, j — 1,2, . . .) in an infinite chain; = \i — j\ is the 
distance between sites and 

This model presents a phase transition at a non-zero critical temperature Tc{a) for 1 < a < 2 [l^ [Tfl| . In the 
ferromagnetic version, Jo > 0, the critical temperature behaves asymptotically as fcfcTc/Jo ~ 2C(a) ~ 2/(1 — a) when 
a ^ 1+ where ({a) = X)^o it^ Riemann zeta function. When 1 < a < 1.5 the critical exponents 

are classical. The thermodynamic limit in the ferromagnetic model is not defined when < a < 1 and the system 
presents non-extensive behaviour^; however, with an appropriated regularization procedure it can be shown that 



3 



the system still presents a phase transition where all the thermodynamic functions (not only the critical exponents) 
obey mean-field behaviourp, E|- These results generalize to ferromagnetic models in arbitrary dimension d with 
a = d + a: the thermodynamic limit is well defined for cr > 0. The critical exponents are classical for < cr < d/2 
and assume continuously varying values for d/2 < cr < 2. When a > min{d, 2) the system exhibits short range critical 
behaviour [l^. 

In this section we considered a finite chain with 1 < a < 2; the finite size version of the Hamiltonian Q is 



H 



5,; = ±1 



(3) 



where now the sum j) runs over all distinct pairs of sites — 1,2, ... ,L) in a chain of lenght L, rij = 

\i — j\ = 0, 1 ... L — 1 and W{rij) is the effective interaction between spins resulting from the different types of 
periodic boundary conditions. 

The FIPBC effective interaction Wp is then simply defined as 



Wpir) = f (min[r, \r — L\ 
while the IPBC effective interaction Wj is defined as 



(4) 



and 



n=0,±l, 
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(6) 



where n — ±1, ±2, . . . , ±00 numerates the infinite repHcas of the system n = 0. Equation lOl is a slowly convergent 
series and cannot be truncated in a straitforward way. In order to obtain a controllable approximation we rewrite 
Eq.lOl as: 
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Since r < L, for a > 1 and m ^ 1 we can approximate Wi{r) « Wmir) for r = 1, . . . , L — 1, where 



(7) 



E 



Ln— ibl....,ibm 

J_ 2C{a) J_ 



1 0^1 
h 2 > — 

In + r/Ll" ^ 



E 



n— ±l,...,±7n 



1 1 



|?i + r/L|" 



(8) 



The finite approximation Wmij) can be calculated once for ever at the beginning of the simulation. In order to control 
the accuracy of the approximation, we can calculate the differences 



AW™(r) EEW^„+i(r)-VF™(r) 
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Noting that 
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a 


IPBC L ~ 100 


Ref.fg] 


1.1 


0.049 ± 0.003 


0.0476168 


1.2 


0.092 ± 0.004 


0.0922314 


1.3 


0.141 ±0.005 


0.136110 


1.4 


0.185 ±0.005 


0.181150 


1.5 


0.238 ± 0.006 


0.229155 



TABLE I: Inverse critical temperature Kc — Jo/Tc obtained with L ~ 100 IPBC compared with the best estimates from 
Ref.[9] 



AWmir) < AW„^iL) - ^ ^ 



L"\(m + 2)" to" (to±1)^ 
we can estimate the error of the approximation by choosing m such that AWm(L) < e; since for to ^ 1 

L°'{m + 1)" 

we have 



?^ ± O (to-) 

(to +1)2 ^ ' 



. a a + 2 
TO > - 1 

In all the simulations we choose to such that e < 10^^ for every value of L and a. 

We first considered the ferromagnetic case Jq > for 1 < a < 2. We calculated the critical temperature using both 
FIPBC and IPBC by means of the so called Binder cumulant|2^ 

C/l = 1 - ^ % 

where M is the order parameter (magnetization per spin). The critical temperatures were estimateded as the point 
of intersection|j22i| of the curves Ul' , Ul" vs. T (in what follows we choose fc^ = 1), for pairs of values {L',L") chosen 
around some average L, with a dispersion around L of less than 10%. The results for boh type of boundary conditions 
were calculated using the same single spin-flip Monte Carlo dynamics (Metropolis). The results are shown in FigQ] 
where the exact asymptotic behaviour when a ^ 1 is also shown for comparison. For a > 1.5 the results are very 
similar, with an obtained difference between both types of boundary conditions of less than 5% and decreasing when 
a ^ 2. On the other hand, for a < 1.5 the difference in performance increases dramatically as a decreases. In table 
ID the results for the best estimate using IPBC for a < 1.5 are compared with the best estimates in the literature 
obtained by Luijten and Blotejp, using a specially designed block Monte Carlo algorithm for ferromagnetic systems 
with long range interactions ][21| (IPBC are implicit in this algorithm) and chains of very large sizes (L = 150000). 
While the IPBC gives a rather accurate result, even for a modest size of L - 100 (see table |l}, the FIPBC with 
L ^ 1000 give an error that reaches around one order of magnitude for a = 1.02 (see FigQJ. 

We next considered the antiferromagnetic case Jq < for 1 < a < 2 and system sizes L = 300, 500 and 1000 using 
both type of boundary conditions, and for < a < 1 using FIPBC (IPBC cannot be implemented in this case). While 
the ferromagnetic version of this model has been extensively studied, to the best of our knowled ge t he only existing 
results for the antiferromagnetic case for 1 < a < 2 have been obtained by Romano (see Ref.[2^ and references 
therein). In particular, no previous results are known for < a < 1, while the a = case has been exactly solved (23|. 
We found that the system presents a low temperature ordered antiferromagnetic phase for every value of < a < 1 
below a flnite critical temperature. At variance with the ferromagnetic case, the antiferromagnetic model shows a 
well defined thermodynamic Hmit for < a < 1; this fact has also been observed in an antiferromagnetic model with 
long range interactions defined in the hypercubic cell|i26y|. 

For this system the critical temperature was estimated as the maximum in the staggered susceptibility, defined as 
the fluctuations in the staggered magnetization. The results are shown in Fig|2| We see that both type of boundary 
conditions give in this case the same result (within the error bars) for every system size when 1 < a < 2. The 
results for L = 1000 are in good agreement with previous estimations [l^l using IPBC for a — 2: Tc/\Jo\ ~ 0.219, and 
converge to Tc = for a ^ 0, reproducing the exact solution(23|. 
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FIG. 1: Critical temperature Tc/Jo vs. a (fcs = 1) for the one-dimensional ferromagnetic model obtained with FIPBC and 
IPBC for different values of L. The notation L ~ indicates the average value of the pairs of values {L',L") used to calculate 
the Binder cumulant. In the inset it can be seen the same figure with a magnified vertical scale. 



Finally, we considered an Ising model with competing short range ferromagnetic interactions and long range anti- 
ferromagnetic interactions defined by the Hamiltonian 

H = -JoY^ S,S, + JdY. /('^y ) ^^^J = ±1' (10) 

<i,3> (ij) 

where Jo, Jd > 0, the sum J2<ij> runs over nearest neighbor paris of spins and /(r) is given by Eq.JSJ. This is 
a one dimensional generalization of a model for an ultrathin magnetic film with competing ferromagnetic exchange 
interactions and dipolar interactions that will be described in the next section. The finite temperature phase diagram 
of this model when 1 < a < 2 is similar to the one corresponding to the two dimensional case. At low temperatures it 
presents different modulated phases (with zero global magnetization) of period 2h, depending on the ratio S = Jo/ Jd, 
with a finite critical temperature Tc = Tc{a, S). The period of the modulation increases with S, ranging from h = 1 
(antiferromagnetic) for low values of S, up to the system size for very large values of S, passing by all the possible 
values commensurable with the system lenght L. The calculation of the global phase diagram of this model is in 
progress and the details will be published elsewhere. 

For the comparison of the boundary conditions we choose some representative values of a and S; other values gave 
similar results. The results are shown in table lUl No significative differences were observed between both type of 
boundary conditions. 



III. TWO-DIMENSIONAL ISING MODELS 

In this section we considered Ising models in a square lattice of sites with interactions that decay as 1/r^. IPBC 
were implemented in this case using the Ewald sums technique^. We started with a Hamiltonian of the type ^ with 
a = 3 and Jo > 0. The choice of the vale a = 3 obeys two reasons: first, it allows the usage of the Ewald sums for 
the implementation of the IPBC and second, because of its close relationship with dipolar interactions. To perform 
the comparison between boundary conditions we calculated in this case the specific heat 
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FIG. 2: Critical temperature Tc/\Jo\ vs. a {ks — 1) for the one-dimensional antiferromagnetic model obtained with FIPBC 
and IPBC for different values of L. 



CiT) = ^ {{H^), {H)l) 

as a function of T for different system sizes L. In FigEI we show the specific heat calculations for both types of 
boundary conditions for system sizes L = 24,32,48,64 and 96. Although FIPBC display larger finite size effects for 
small system sizes than the IPBC, both type of boundary conditions give similar results when L = 96, as can be 
appreciated in Fig^l where we compare the results of FIPBC and IPBC for the smallest and the largest sytem sizes. 
For L = 96 the curves obtained wiht both kind of boundary conditions show clearly the emergency of a discontinuous 
shape, consistently with the expected classical behaviourigj. 

Finally we considered the two-dimensional version of Hamiltonian ifTTill . This model represents an ultrathin magnetic 
film where the axis of easy magnetization is oriented perpendicular to the plane of the filmQ. In this case the first term 
in Hamiltonian represents ferromagnetic exchange interactions while the second represents dipolar interactions. 
The system presents a low temperature ordered phase below some critical temperature Tc = Tc{5), composed by 
ferromagnetic stripes of width h = h{5), so that spins belonging to adyacent stripes are antialigned. The overall 
known features of the finite temperature phase diagram of this model are described in Refs.0-0|. For the comparison 
of the boundary conditions we choose the value (5 = 3, for which the low temperature phase corresponds to an ft. = 4 
striped phase, since it is representative of the global behaviour for intermediate values of 5 and was previously anaHzed 
in detail by Booth et alQ using Monte Carlo with IPBCjU. The results are shown in FigElfor L = 24,32. We see 
that both types of boundary conditions give almost the same results even for rather small sizes, except very near the 
critical point (left peak; see Ref.|3| for details) where the FIPBC give a little larger peak than the IPBC. 



IV. DISCUSSION 



We have compared the finite size effects introduced by FIPBC and IPBC in one and two dimensional Ising models 
with different types of long range interactions. This allowed us to anahze a large variety of situations of interest and 
the infiuence of different parameters. We have shown that, when the interactions are such that the system presents 
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FIG. 3: Specific heat curves for the two-dimensional ferromagnetic model with a = 3 obtained with (a) FIPBC and (b) IPBC 
for different values of L. 




FIG. 4: Comparison between the specific heat curves for the two-dimensional ferromagnetic model with a = 3 obtained with 
FIPBC and IPBC. 
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1.2 


0.076 ± 0.002 


0.075 ± 0.002 
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0.060 ± 0.002 


0.0625 ± 0.002 


0.5 


1.6 


0.0575 ± 0.002 


0.060 ± 0.002 
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0.0625 ± 0.002 


0.0675 ± 0.002 
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0.0575 ± 0.002 


0.058 ± 0.002 
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0.0625 ± 0.002 
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0.061 ± 0.002 
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0.065 ± 0.002 


0.065 ± 0.002 
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0.077 ± 0.002 


0.075 ± 0.002 


1.3 


1.4 


0.073 ± 0.002 


0.073 ± 0.002 


1.3 


1.6 


0.074 ± 0.002 


0.073 ± 0.002 


1.3 


1.8 


0.0725 ± 0.002 


0.073 ± 0.002 



TABLE II: Critical temperature T^/ Jd for the one-dimensional model with competitive interactions l|10|l obtained with FIPBC 
and IPBC for different values of a and 5. 
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T/J, 

FIG. 5: Specific heat vs. T/Jd in the two dimensional model with competing interactions l|in|l with (5 = Jo/Jd = 3. 

a low temperature ordered phase with zero global magnetization the finite size effects produced by both types of 
boundary conditions are almost the same. This can be understood if we consider that, in this case, the infinite sums 
in the local fields calculations converge fast and hence the finite size effects associated with the long range character 
of the interactions are wiped out quickly; hence, the truncation associated with the FIPBC is enough to take it into 
account. A similar effect was observed by Kretschmer and Binder[18j in the antiferromagnetic phase of the three 
dimensional Ising model with nearest neighbor exchange interactions and long range dipolar interactions. On the 
other hand, when the system presents low temperature ferromagnetic ordering, large differences in the finite size 
effects produced by the different periodic boundary conditions can be expected, depending on the values of d and a. 
Moreover, our results suggest that the large differences appear when the ratio a/d is such that the system enters into 
the classical regime a/d < 3/2 (cr < d/2), while for a/d > 3/2 both type of boundary conditions produce similar 
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finite size eflects, with a decreasing difference as a/d increases. This is also consistent with Kretschmer and Binder 
result for c? = 3, a = 3. 

We have also shown that in the ferromagnetic one dimensional model the finite size effects given by the FIPBC 
becomes very large when a ^ l"*". This suggests that the same problem could be observed in the non-extensive regime 
< a < 1, where the FIPBC are the most direct form of implementing periodic boundary conditions. Moreover, in 
this regime when a < d the implementation of the IPBC is more involved because the infinite sums of the type © 
do not converge. Since in the a = case we lost the boundary completely (this corresponds, with an appropriated 
re-scaling, to the Curie- Weiss model) we can expect that in the limit a ^ the infiuence of periodic boundary 
conditions will be also negligible. Hence, we may expect the largest finite size effects when using FIPBC in the region 
a/d^ 1. 

Finally, as a by-product of this work, we introduced a simple method for implementing the IPBC in one dimensional 
models for arbitrary values of a > 1. 
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